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I. INTRODUCTION 


The first free electron laser (FEL) was demonstrated in 1976 by Madey and his 
co-workers at Stanford University [Ref. 1]. Since this time, scientists from around the 
world have worked to develop the field of FELs both experimentally and theoretically. 

FELs have become an interest of the United States Navy as a possible solution 
to an important current problem. One of the most real and dangerous threats to Naval 
ships is anti-ship missiles (ASM) which are included in the arsenals of most every 
country in the world. With the increasing mission of single-ship operations and 
operations close to land masses world wide, each US Navy ship needs to be able to 
defend itself against a sudden ASM attack. The current close-in-weapon-system 
installed aboard ships is the Phalanx gun. This weapon is inadequate for protecting 
the ship against incoming missiles because the missile approaches too close to the 
ship before being destroyed by the gun. As a result, the Navy is exploring different 
options for a solution to this problem. 

The advantage of lasers is the speed of propagation of the destructive energy, 
the speed of light, and the ability to focus intense energy to a small spot at 5 km to 10 
km. These characteristics can increase the distance at which incoming missiles are 
destroyed, lessening the likelyhood of high energy debris impacting the ship. The 
flexibility of FELs has resulted in Navy funding to develop a MW FEL. Missiles require 
about 20 kW/cm^ over 100 cm^ to crack the structure which causes a rapid self 
destruction of the high speed missile. Therefore, a MW laser seems to be the range 
of power output suitable for this endeavor. 

Currently, a 1 KW FEL is being constructed at Thomas Jefferson National 
Accelerator Facility (TJNAF) in Newport News, VA. This FEL will be the most powerful 
ever constructed. For this reason, the Navy’s directed energy office is assisting in it’s 
development to gain knowledge on high power FELs and their feasibility for a ship 
self-defense weapon system. 
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In hopes of developing a feasible design, the Navy assembled scientists from 
around the country for a 3 day MW FEL design meeting at TJNAF in September 1996. 
From the many designs presented, two FEL weapon designs emerged. These two 
designs, the MW regenerative amplifier design and the MW oscillator design, are 
reviewed and their strengths and weaknesses discussed. 

Some analysis of each design determines the possible feasibility of these 
proposed FEL weapons. Analysis through computer simulations is conducted to 
evaluate the current parameters of the MW regenerative amplifier. A potential problem 
with the MW oscillator design stems from coherent synchrotron radiation (CSR) and its 
effect on emittance of the electron beam. This problem is reviewed along with the 
conclusions of an electron beam transport workshop held at NFS. Possible CSR 
experiments and parameters are listed and a relative comparison is made. One 
possible experiment involving the Princeton FEL (CIRFEL) is discussed and some 
computer analysis is conducted. This work is intended to aid the Navy’s Directed 
Energy Office in its quest for a ship self defense FEL weapon. 
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II. BACKGROUND 


A. THE BASIC FREE ELECTRON LASER 

An FEL uses a relativistic beam of electrons to generate electromagnetic 
radiation. A transverse, periodic magnetic field, called an "undulator", wiggles the 
electrons causing acceleration when the electrons experience Lorentz forces in the 
undulator as shown in Figure 1. 



Figure 1. This is a typical undulator. The magnetic field 
causes the electrons to wiggle in the transverse plane. 
From Ref. [1]. 


The result is an interaction between the electron beam and the co-propagating 
electromagnetic wave and is the basic operating principle of all FELs. [Ref. 1] 

The basic components of an FEL are an electron accelerator and an undulator. 
There are two distinct FEL designs: the oscillator and the amplifier. The oscillator has 
an optical cavity enclosed by mirrors with the undulator at the mid-section as shown in 
Figure 2. The oscillator transfers a small percentage of energy from the electron beam 
to the light wave on each successive pass of the oscillating light wave in the undulator. 
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electron beam 



resonator mirror 


Figure 2. This is a typical FEL oscillator. The electron 
beam passes through the cavity and is acted on by the 
magnetic field from the undulator. Light is represented by 
the shaded portion between the two curved mirrors. 
From Ref. [2]. 


The amplifier looks similar but lacks the mirrors. The amplifier relies on a high 
percentage of the electron beam energy being transferred to the light wave on a single 
pass through the undulator. Typically, the oscillator will be designed to have low 
single-pass gain and the amplifier will have high gain. 

Many parameters common to all free electron lasers have typical values. The 
distance from the center of one magnet to the center of the next magnet is defined as 
the undulator wavelength, Xq, which varies from 2 cm to 10 cm. The number of 
periods, N, in undulators varies from 20 to several hundred. Since the length of an 
undulator is L = NXq, undulators range in length from 1 m up to 20 m. Electron 
beams used in FELs have radii, r^,, from a fraction of a mm up to a cm. Peak currents 
vary widely in FELs. Some peak currents are around 1 A and some are as high as a 
kA. The average current is much lower and has nominal values in the nA to mA 
range. One of the most desirable attributes of FELs is that they can be designed to 


4 





provide light at specific wavelengths, X. The optical wavelengths of FELs range from 
X, = 1 cm to X = 200 nm. 

B. FREE ELECTRON LASER THEORY 

1. The Pendulum Equation 

The electrons traveling through the undulator experience forces from both the 
optical and undulator fields. By examining the electron dynamics classically, it can be 
shown that the equation which governs the electron’s microscopic motion is the well- 
known equation describing the motion of a classical pendulum. This surprising parallel 
is extremely helpful in the understanding of free electron lasers. What follows is a 
derivation of FEL microscopic electron dynamics. [Ref. 3] 

An important relation in FEL dynamics is the resonance condition. Light 
traveling at speed c overtakes electrons traveling at speed in the undulator, 
where and is the longitudinal speed of the electron. The definition of 

resonance is when one wavelength of light X passes over an electron as the electron 
travels through one undulator wavelength Xq . The time it takes an electron to travel 
through one undulator period is t = Xq/v^ = Xo/(P 2 c). After one undulator period, the 
one wavelength of light will be ahead of the electron a distance X, where 
X= ct - v^t = cf(1-P2). Substituting for time gives the resonance condition, 

X = Xo^^-^ . (1) 

Pz 

For highly relativistic electrons, as in an FEL, P^ = 1- In this case, the wavelength of 
the light will be much less than an undulator wavelength, X «; Xq. 

A static helical field from the undulator magnets takes the form 

= B{tcosl<oZ + ysin/foz) , (2) 
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where /Cq = 2hJXq, and z is the longitudinal position along the undulator. The undulator 
wave number is /cq and Xq is the undulator wavelength. B is the amplitude of the 
undulator magnetic field and /, j, and /?, are unit vectors in the x, y, and z directions 
respectively, as represented in Figure 1. (Note that the field in Figure 1 represents a 
linear undulator, not one that is helical.) The magnetic field strength for undulators is 
typically 2 to 7 kG. A circularly-polarized plane-wave is assumed for the optical field, 

= E(fcos\y - ysin\y) , = E(/sin\j/+/cosxy) , (3) 


where \|/ = /cz - tof + <(>, and k = co/c. The electric and magnetic field amplitudes are 
E, k is the optical wave number, co is the angular frequency of the optical wave, c is 
the speed of light in vacuum, t is time, and ^ is the optical phase. 

The Lorentz force equations are. 


dt 




djymc^) 

dt 


= -e V 


(4) 


where the electron velocity is V = ^c, and the Lorentz factor is y = (1 - The 

electron momentum is p = ymV, the electron energy is ymc^, m is the mass of an 
electron, and e is the electron charge magnitude [Ref. 4]. The electron energy in 
FELs ranges from a few MeV with y-S all the way to a few GeV with y-5000. 
Substituting for the momentum in the Lorentz equations gives, 

dt me l 



where ^ = iu + • 

The electric and magnetic fields, (2) and (3), are then inserted into (5) to yield. 


dt 

djy^z) 

dt 


= —— E(1-p,)(fcos\y - ysinxy) + BP 2 (-rsin/foZ + ycoskoz) , (6) 

me L J 


e 

me 


E(P;(.cos\|/ - PySin\|/) + e(PxSin/foZ - PyCoskoZ) 

-^ =--^E p;,cos\K - pysinvj , 


l< , 


(7) 

( 8 ) 
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where = (3;^^ + Py/, and - Pz^- It is assumed that for relativistic electrons where 
p^ = 1 , E (1 - P 2 ) c 6 p^, so the first term on the right-hand side of ( 6 ) is ignored and 
P^ is set equal to one in the second term. 

Remembering that z = ct very nearly and solving (6) for gives 
P^ = -{K/i)coskQZ and Py = -{K/y)s\nkQZ, where the undulator parameter is 
K = eB/{kQmc^). From the typical FEL parameters, K=^. Perfect injection into 
helical orbits is assumed so the constant of integration is zero. Substituting the result 
for into (8) gives the following equation for the electrons, 

-*=T=-^cosK + « , (9) 

ut ymc 

where the electron phase in the combined undulator and light fields is 
= {k + /fo)z - cof, and d/dt = ( ) . Squaring the result for gives pj = 

which further yields 1 / 7 ^ = ! - K^/ 7 ^ - p| from the definition of the Lorentz factor. 
Solving for p^ and substituting into ( 1 ) for highly relativistic electrons where y» 1, the 
resonance condition becomes. 




Xo(i + k^) 

2f 


( 10 ) 


Now the resonance condition is in terms of the electron energy. [Ref. 5] 

To examine take time at f = 0 so ^(0) = {k+ko)zo = {2nlX)Zo, where Zq is the 
initial electron position. We know Jlq which also means that k » /Tq- Changing 
the electron initial position by an amount ranging from zero to one wavelength of light 
causes ^ to vary between 0 and 2jt. From (9), a change in ^ of 71 causes y, the rate of 
change in energy of the electrons, to reverse sign. Since the electron beam spans 
many wavelengths of light, the electrons are uniformly distributed throughout each 
wavelength. For this reason about half the electrons will gain energy and about half 
will lose energy. As electrons gain energy, their velocity increases while electrons that 
lose energy decrease in velocity. This leads to bunching of electrons. Since y > 0 for 
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about half the electrons and y < 0 for the rest, there does not seem to be an overall 
change in energy, and therefore no gain. Gain will be discussed in the next section. 

The next step in finding the pendulum equation is to relate y to using 

y"^ = 1 - p^ = 1 - P| - p|. The values of p^ and py found previously give 
(1 + K^)y“^ = 1 - p| ■ Taking the time derivative of each side and rearranging gives. 


y Y^PzPz 
Y (1 + K^) 


( 11 ) 


The final step in arriving at the pendulum equation is to relate y obtained in 
(11) to Using the definition for the second derivative with respect to time is 
^ = (/f + I<q)c^z . This can be substituted into (11) to yield a relation between y and 
Since free electron lasers use highly relativistic electrons, p^ = 1, and /f » /Cq so that 
the relation becomes y/y = y^^/[(ji)(1 + K^)]. Substituting the resonance condition (10) 
into the newest relation between y and ^ yields y/y = ^/2(0o . Finally, substituting in the 
value of y in (9) gives the microscopic electron equation of motion as. 




2(iiQeKE 

■f-mc 


cos(C + ({)) . 


( 12 ) 


Now the microscopic electron phase dynamics are in the form of the pendulum 
equation. [Ref. 3] 

Dimensionless parameters are introduced to enhance the understanding of 
FELs. The first such parameter will represent time. For an average electron, the time 
it spends in the undulator is determined by the length of the undulator L and the 
speed of the electron Vq such that t = L/Vq . Remembering that Pq = Vq/c then 
t = L/PoC = L/C since, for highly relativistic electrons, po = 1 . Now the dimensionless 
time T will be defined as t = ct/L . 


Substituting this dimensionless time parameter into (12) yields, 

00 

; = lalcos(C + (t>) , (13) 

o 

where la I = 47t/Ve/CLE/(y^mc^), the number of undulator periods is N = L/Xq, ( ) 
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represents differentiation with respect to x, the dimensionless optical field is 
a = lale"**, and / = Typical values of the dimensionless optical field amplitude 
are la I for weak optical fields and la I > Jr for strong optical fields. The 

assumptions that have been made are y» 1 for highly relativistic electrons so that 
^0 » X, and therefore k » /Cq. Examination of (13) shows that for -k/2 < (^ + ()>) < k/2, 

oo 

the phase acceleration ^ is greater than zero, and for 7t/2 < (^ + <j)) < 3k/2, it is less 
than zero. Electrons that begin their evolution along the undulator randomly spread in 
phase space from -ti/2 to 3 ji/ 2, begin to bunch at (C + ())) = k/2 because of the 
differences in phase acceleration between the electrons surrounding ti/2. 

As the electron beam and optical wave interaction evolves, the optical field 
strength passes from the weak field regime to strong fields. In weak optical fields 
la I «: 71 , bunching of electrons is very slight. In strong fields la I ^ ti, electron 
bunching occurs rapidly and the electrons begin to become trapped in closed phase- 
space orbits, the boundary of which is represented by the "separatrix." The bunched 
electrons that are trapped, after causing gain, then shift in phase to one which causes 
the electrons to begin to absorb energy from the optical field, the result of which is a 
decrease in the optical field gain. As the gain begins to diminish the laser is said to 
have reached saturation. Gain is explained in detail in section 3 of this chapter. 

An example is seen in Figure 3 where the final positions of 3000 sample 

o 

electrons are shown in phase space. The electron phase velocity v, defined as v = ^, 
is plotted on the vertical axis, and electron phase ^ is plotted on the horizontal axis. 
The electrons’ initial phase velocities are distributed in a Gaussian with standard 
deviation ©6 = 1 around Vq = 3. A number of electrons are trapped in closed phase- 
space orbits inside the solid line representing the separatrix. This is typical of a low 
gain oscillator in strong optical fields near saturation. The gain is plotted at the upper- 
right on a logarithmic scale from the beginning to the end of the undulator represented 
by X = 0 to X = 1, respectively. Note that saturation has been reached, marked by the 
decrease in gain near x = 1. Strong optical fields are represented since the initial 
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optical field strength ao = 20> ti. The dimensionless current density y = 10 and the 
dimensionless energy spread Cq = 1 are explained further in section 3 and 4, 



Figure 3. A plot of FEL phase space evolution for low 
gain. 


respectively. There are N = 25 undulator periods and the optical phase is also shown 
in the lower-right portion of the figure. The optical phase is initially zero and steadily 
increases up to a final value of <[) = 0.15. 

2. The Optical Wave Equation 

The pendulum equation, (13), governs the motion of the electrons in the 
undulator. The optical wave equation will now be developed so that the evolution of 
the complex radiation field can be described mathematically. The wave equation in a 
vacuum with no current density is of the form 
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A{x,y,z,t) = 0 


(14) 


’^2 _ JL Ji' 

c 2 dt\ 

where ^ = f d/dx + j d/dy + k didz and X is the vector potential from which ^ and ^ 
may be derived by ^ = (- 1 /c) d^/dt, and ^ = ^x>^. A vector potential is selected in 
the form of a circularly polarized plane wave, 

X = [/sin\|/ + ycos\|r] 

K 

where E is the complex electric field. The vector potential can be written in the 
following form 


k 


(15) 


where e= IE( 5 ?,f)le'‘*’^^’^\ ^ represents the optical phase, and 6 = (-/, 1 , 0 ) is the 
polarization vector. When evaluating (14) by substituting in (15), the slowly-varying 
amplitude and phase approximation is made. In this approximation, the field is taken 
to vary slowly in time during an optical period (E «: coE , (j) «: ox))), and slowly in space 
over an optical wavelength (E c kE , k<^) where (') represents differentiation 

with respect to z. This is an appropriate assumption since an FEL emits a relatively 
narrow spectrum or linewidth. If the fields varied quickly, the emission spectrum would 
be broadband which is not the case. [Ref. 5] 


Using e(x,y, 2 ,f) and evaluating (14) by using (15), yields 

8 (x,y,z ,0 = 0 


+ ik 

A. + IAI 

p 4- 

dz C dt 


(16) 


where terms of e and e" have been ignored because of the slowly-varying amplitude 
and phase approximation, and ^ 4 ^= t d/dx + j d/dy. The first term in (16) describes 
diffraction. A solution to the second term alone is of the form e= f{z-ct) where f 
represents any function. Using t = ct/L and defining 2 = z-ct, it can be found that 
d/dz -I- (1/c) d/dt = (1/L) d/dx where the 2 terms have cancelled. In dimensionless 
parameters the optical field is represented by a = lale'**’ which is of the same form as 
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e=IEle'‘'’. The dimensionless optical field amplitude la I = 4nA/e/<L IEI/(7^/77c^) 
includes the electric field amplitude. Exchanging e for a, and including t gives 


.72 + 


L Bx 


a =0 


(17) 


Again, the first term describes diffraction. If the two terms are comparable, then 
diffraction is important. If the first term is small in comparison with the second 
because the wave front is very wide, diffraction is not as important and can be 
ignored. 

To obtain a completely dimensionless wave equation, two dimensionless terms 
are defined as ;<: = xV/f/(2L) and y = yV/f/(2L). Using these terms = (/f/2L)Vf 
where + d^ldy^. Using this in (17) gives the dimensionless parabolic wave 

equation 


—V|+ ^ 

4 dx 


a = 0 


(18) 


This derivation began assuming no current density was present. In an PEL with an 
electron beam (18) becomes 




a = -j<e 


(19) 


where = {k + kQ)z - (at is the electron phase in the combined undulator and light 
fields, and j is the dimensionless current density, j = 8N{e%KL)^p/{’fmc^) described 
in detail in the next section on gain [Ref. 5]. Values of J vary a great deal from one 
PEL to another with an approximate range from unity to several thousand. 

In an oscillator PEL, the curvature of the resonator mirrors and diffraction play 
important roles in determining the shape of the optical mode. Spherical mirrors with 
y = 0 give rise to Gaussian optical modes and the solution to (19) is called a 
"Gaussian" beam because of the optical mode’s transverse shape. This solution with 
y = 0 is 
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( 20 ) 



and 


(t)G(x) = -tan ^ 


Zo + (t - 


^ - '^w) 


where w^(t) = 1 + (x - t^^)^/Zo , r = (x^ + is the radial coordinate, Zq is the 
dimensionless Rayleigh length, <})g(x) is the phase shift of the Gaussian mode, and x^^ 
is the position in x of the optical mode waist. For example, if Xj^ = 0.5 then the optical 
mode waist is located at the midpoint of the undulator. The Rayleigh length is given 
by kWq/X where Wq is the size of the optical mode waist. The Rayleigh length is the 
distance in which the initial optical mode area doubles in size. The dimensionless 
Rayleigh length is then defined as Zq = %WqI{L'K). This Gaussian beam represents an 
optical mode which contracts in the middle and expands at the ends and depends on 
the spacing and curvature of the mirrors. If the optical mode is a Gaussian, the initial 
phase velocity that yields resonance is shifted to Vq = I/Zq instead of Vq = 0 when Zq is 
large. [Ref. 6] 

In some cases, diffraction may not be important so that (19) simplifies to 


O 

a = -i<e ‘^> 


( 21 ) 


Now that the equations which govern the electron dynamics (13) and the optical wave 
(19) have been developed, the evolution of how the optical fields grow will be 
discussed. 

3. The Gain Equations 

FELs can be designed to have either low gain, as in a typical oscillator design, 
or high gain, as in a typical amplifier design. Gain is simply the fractional change in 
power of the optical field in a single pass through the undulator. From the 
conservation of energy, an energy change in the electron beam must result in an 
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opposite energy change in the optical field. The gain may be calculated using the 
change in energy of the electron beam divided by the optical wave energy. 

The first step in in finding the average change in energy of the electron beam 
will be to find the change in energy of a single electron, Aymc^. This change in 
energy will depend on the change in phase velocity of the electrons with respect to 

o 

time. The phase velocity is defined as v = ^ and is v = L[(/f+/fo)Pz-^]- Now a relation 
needs to be found between Ay and Av. Remembering that k » /Cq. then Av = LkA^^. 
Now using N = L/Xq, X = 2n/k, Xq = Znlk^, and (10) gives Av = N2rt{2'^/{1+K^))A^z- 
Using (11), Av gives. 


Av = 4kN-^ . 

Y 

The change in energy of an electron becomes Aymc^ = ymc^Av/{4TiN). 


( 22 ) 


The average change in energy of the electron beam, the change in phase 
velocity, Av, is the difference between the average phase velocity, <v>, and the initial 

o 

phase velocity, Vq, of the electrons, Av = <v> - Vq,. The relation v = ^ will be used to 
eventually obtain <v>. Using (13) in weak fields, ^ can be expanded in powers of the 
optical field ao so that ^ + • • • . For the expansion, a(0) = a© and 

v^°^ = Vq. So now = Vq and integrating with respect to t will give = vqX + ^o- 

oo 

Now (13) will become ^ = aocos(^^°^) = aocos(^o + Successive integrations will 

yield so that 


C = Co + voi: - ^ [cos(Co + vot) - cos(Co) + VoT:sin(Co)] + • • • • (23) 

All the terms in the expansion of C can be found by repeating this sequence. To 
obtain v, only terms of and terms of to first order will be kept ignoring (C^^^)^ 
and higher order terms. Differentiating C gives 


14 


(24) 


V = + v<2) + = Vo + — 

Vo 


sin(;o + VoX) - sin(;o) 


^2 


J_ 

4 


cos(2Co + 2vox) - cos(2Co) + cos(vox) - 1 - VoXSin(Co)cos(Co + Vox) 


+ . . . 


The next step will be to obtain <v> from (24). 

The average changes in phase velocity of an electron will be obtained by 
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(1/2tc)|vc/^o • Using this to evaluate (24), the last term in (24) will be the only term to 
survive so that 

<v> = + • • ■ 


r 

"0 


2Vn3 


2cos(vox) -2 + voxsin(vox) 


+ 


(25) 


This relation will be used to find Av in (22). Now the average change in energy of an 
electron is 


Aymc^ = 


Ymc^(<v>-Vo) 
471 A/ 


(26) 


Before introducing the gain equation for low gain, the "filling factor" F, will be 
introduced which is the area of the electron beam divided by the area of the optical 
beam, F = ! Wq, where r^, is the radius of the electron beam and Wq is the optical 

beam waist. This factor arises because only a fraction of the optical mode is filled with 
electrons and accounts for the fraction of the two beams which interact with each other 
and therefore affect the gain. Now the average change in energy of the entire electron 
beam will be (pFdV’)(Aymc^) where p is the electron density and dV is the volume 
occupied by the optical mode one wavelength of light long. A typical value of 
p = 10^^ cm"^ for a current of / = 100 A and r^, = 1 mm. [Ref. 7] 

The radiation energy contained in a volume dV is 2E^dV/8K [Ref. 7], so the 
gain equation is 

Q [pFdV][ymc^{<v>-Vo)/{4KN)] 

2E^dV/8K 
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This equation is true only for low gain since Sq has been held fixed when interpreting 
(23) and (24). In weak optical fields, ^ can be expanded to obtain <v>, substituting 
(25) into (27) and simplifying yields the following equation for low gain in weak fields, 


G = j 


[ 2 - 2cos(vox) - voxsin(vox)] 


(28) 


In (28), j = 8N{eKKLfpF/{-^mc^) is the dimensionless current density, and is the 
most important FEL parameter. Low gain corresponds to y ^ and high gain 
corresponds to / » tc . 

An example of low gain in weak fields is shown in Figure 4. The figure shows 
the final gain spectrum G(vo) at x = 1, the end of the undulator, plotted against the 


**** Gain Curve **** 

j=l a =0.01 N=25 

o 

Maximum gain at v^= 2.6 

0.13 

0.0 

-0.13 

-12 V 12 

o 

Figure 4. The gain spectrum G{vq) for low current and 
weak fields. 



initial phase velocity Vq. The dimensionless current density y = 1 < n is in the low gain 
regime and the initial optical field strength Sq = 0.01 «: k represents a weak field. 
There are N = 25 undulator periods with a maximum gain of G=0.13 at phase 
velocity Vq = 2.6. The gain is anti-symmetric about the resonance Vq = 0, and positive 


16 










gain occurs for vq > 0. For vq < 0 the gain is negative which means the electron beam 
is absorbing energy from the optical wave. A shift in vq = 2k can cause a shift from 
positive gain, amplification, to negative gain, absorption. 

The development of the gain equation for large /' » ji requires the use of the 
wave equation, because the optical field grows rapidly. The expansion performed to 
obtain the low gain equation cannot be used. For high gain, the governing equation is 

G(,) = leW'2)">'5.l , ,29) 

from [Ref. 8]. The gain is exponential in x along the undulator with growth rate 
proportional to For high gain, the growth rate can be very high. 

An example of high gain in weak fields is shown in Figure 5. The figure shows 
the final gain spectrum G{vq) at x= 1, the end of the undulator, plotted against the 
initial phase velocity vq. The dimensionless current density / = 1000 » it is in the 
high-gain regime with N = 200 undulator periods which is typical of an FEL amplifier. 


★ ★★★ Gain Cuzi/’e ♦*** 

j=1000 a =0.01 N=200 
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Maximma gain at v^= 0.5 

4.2x10^ 


0.0 

-16 v„ 16 

o 

Figure 5. The gain spectrum G(vo) for high current and 
weak fields. 



17 






The maximum gain of G= 4.2x10'* occurs near resonance, Vq = 0 and the gain 
spectrum is almost symmetric about resonance. 


4. Electron Beam Quality 

One of the basic components of every FEL is some type of electron gun and 
accelerator. This initial production of electrons is important to the resulting output of 
the FEL. There is a tradeoff behween high beam current and good beam quality: 
increasing beam current often comes at the expense of beam quality. The 
dimensionless current /oc/A/^ where / is the actual current and N is the number of 
undulator periods, so increasing the number of undulator periods would be an 
approach to increasing j [Ref. 9]. However, an FEL with a longer undulator is more 
sensitive to beam quality than a shorter undulator because of the narrow gain 
spectrum bandwidth. 

To represent the angular distribution of electrons about the z-axis, Oe is defined 
as 


O0 = 


2nN'f^ 

1 + K'2 


(30) 


where 6 is the root mean square (rms) injection angle of an electron. The distribution 
of the electrons in position is represented by 


Oe 



(31) 


where Tq is the rms distance of the electron from the z-axis. [Ref. 9] 

A term commonly used to describe beam quality is emittance. The rms 
emittance is defined as e^ms = ^ 6 Tq. It turns out that the emittance is a constant 
along the length of the undulator [Ref. 10], so decreasing the electrons injection angle 
0 results in an increased displacement from the z-axis, Tq. Using (30), (31), and (10) 
the definition of emittance becomes e^ms = where 
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is the dimensionless emittance. Another useful definition of emittance is normalized 
emittance defined as 


~ y^rms > 


(33) 


which is usually stated in units of n mm-mrad. 

While emittance determines beam quality in the transverse direction, it does not 
correspond to the beam’s longitudinal distribution. The term which describes this 
longitudinal spread is 


Oq = AnN-^ 
Y 


(34) 


where Aymc^ is the rms change in energy of an electron from the resonance energy 
[Ref. 9]. The longitudinal spread corresponds to an energy spread in the electrons. 

Electrons do not all enter the undulator perfectly but are randomly distributed in 
both the longitudinal and transverse directions. A Gaussian distribution of the 
electrons is assumed of the form 



(35) 


where a is the standard deviation and x is the parameter which is being distributed 
[Ref. 9]. An example of the effect of beam quality on an PEL is shown in Figure 6. As 
with Figure 3, this phase space evolution is generated by solving the pendulum 
equation (13) and the wave equation (21) numerically. The dimensionless current 
density y = 10 is at the low end of the high gain regime and strong fields are 
represented by the initial optical field strength ao = 20 > n. The initial phase velocity 
is vq = 3 with N = 25 undulator periods. The longitudinal energy spread Oq =4 is 
increased from Oq = 1 in Figure 3. This degradation in beam quality causes less 




bunching which is visible when comparing Figure 3 to Figure 6. The energy spread is 
responsible for a decrease in final gain to G = 0.1 and a smaller optical phase shift ((). 


*★* pel Phase Space Evolution *** 

j=10 a =20 V =3 N=25 

o o 
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Figure 6. A plot of FEL phase space evolution for low 
gain with electron beam degradation. 


While gain is the fractional increase in the optical field la 1^, the FEL efficiency 
is the fraction of the electron beam energy extracted by the optical wave. The 
"natural" efficiency of an FEL as it reaches saturation is ti * = 1/(2A/). If N is increased 
to increase gain, the efficiency is decreased. Ultimately, the true efficiency of an FEL 
is simply the ratio of the power taken from the electron beam compared to the initial 
power in the electron beam. [Ref. 11] 
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5. The Tapered Undulator 

As an FEL reaches saturation the electrons have given up enough energy so 
that they are no longer in resonance with the optical wave and gain decreases. 
Tapering an undulator changes the undulators’ resonance condition along the length of 
the undulator so the electrons remain in resonance longer, hence increasing the 
overall gain in strong fields. Tapering can be done by decreasing the undulator 
wavelength Xq in the z-direction, or decreasing the magnetic field strength B, or both. 
Tapering essentially has the same result as accelerating the electrons. So, tapering 
acts as an artificial acceleration to the electrons. This artificial acceleration is 
represented by a dimensionless parameter 6 . When Xq is decreased, this artificial 
acceleration is given by 5 = -ZizN^Xq/Xq and when the undulator field strength is 
decreased it is given by 5 = -4KNK^AB/B{^+K^). Some criteria for tapering must be 
met and are summed up in the following relations as 

lal >6> 4laP'2> 2 k . (36) 

The tapered undulator is more efficient in strong optical fields than an untappered 
undulator. This is where the taper is designed to be effective. However, a tapered 
undulator has lower gain in weak optical fields than an untapered undulator. The 
efficiency of a tapered undulator is estimated by as TI 5 = 5/{8kN). For example, if 
5 = 20jr and A/= 100 periods, then TI 5 = 5/(2A/) while the natural efficiency is only 
Ti' = 1 /( 2 A/). So tapering, in this case, has increased the efficiency by a factor of 5. 
Tapering is a way to obtain higher overall gain in strong fields and higher efficiency. 
[Ref. 12 ] 

An example of the tapered undulator in strong optical fields is shown in Figure 
7. This example is based on a regenerative amplifier design developed by the FEL 
group at NFS and presented at the MW workshop September 24-26, 1996, in Newport 
News, VA. A regenerative amplifier FEL is basically the same as an amplifier FEL. 
The difference is that a regenerative amplifier takes approximately 1 % to 10 % of the 
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optical field (0.01% to 1.0% of the optical power) from the end of the undulator and 
feeds this back to the beginning of the undulator. The light that is fed back is 
focussed at the front of the undulator. This design is based on an undulator length 
/. = 8 m, peak current of 600 A, and a goal for gain G ^ 1000. The dimensionless 


*** FEL Phase Space Evolution *** 

j=6000 a =60 V =0 N=200 

o o 

0^=2 5=160jr 
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Figure 7. A plot of FEL phase space evolution for high 
gain in strong fields with a tapered undulator. 


current density j = 6000»Ji is in the high gain regime and strong fields are 
represented by Sq = 60 > n. As seen in Figure 5, the phase velocity for maximum 
gain occurs at resonance, Vq = 0. The number of undulator periods, N = 200 is typical 
of an FEL amplifier and the Gaussian spread Og = 2 is nominal. The tapering 
parameter 5 = 160;: turns on at x = 0.3 and results in final gain of G = 1000 and 
extraction efficiency ti = 11%. Gain is roughly exponential from x = 0.1 to x = 0.3 
where the taper turns on, and is represented by a straight line on the log plot. 
Saturation is then reached, but tapering allows gain to still increase slowly instead of 
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decreasing, after saturation, as in Figure 3. The trapped electrons are responsible for 
gain and are shown inside the separatrix. The untrapped electrons near v = 400 are 
shifted in phase velocity by the tapering whereas the trapped electrons remain near 
resonance because of strong fields. The optical phase shift Acj) = 1 .Srr over the length 
of the undulator is characteristic of high gain. This phase shift is plotted in the lower 
right-hand portion of the figure. A value of <]) > rc is represented by a negative value of 
phase. The actual value for <[) in this example simply continues to grow from (j) = rt to 
<|) = 1 .5k. 

Normally, natural diffraction of light would cause the optical field to expand in 
the transverse direction over the long propagation length of this undulator. If this 
occurred, the electron beam would become uncoupled from the optical field and the 
FEL interaction would stop. High current density of the electron beam, as seen here, 
causes light to focus back down on the electron beam, hence counteracting diffraction. 
The phase shift imposed by diffraction is opposite to the phase shift resulting from the 
high current. This effect is termed optical guiding [Ref. 13]. The large optical phase 
shift in Figure 7 represents optical guiding which is an important part of a high gain, 
high current FEL. 
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III. MW-CLASS FREE ELECTRON LASER 


A. FEL PARAMETERS 

In order to destroy an incoming missile, the laser power required on the missile 
surface is about 2 kW/cm^ over a ivq = 6 cm radius spot for a 3 second duration. The 
extinction coefficient due to aerosols at sea level is a = 0.05/km, and e~^ estimates 
the power remaining at X = 1 pm wavelength over a distance z. For an example, 
assume an oscillator FEL with N = 25 periods, a peak current J = 600 A, a pulse 
duration of 3 ps, and a repetition frequency of 500 MHz. For X = 1 pm with a nominal 
undulator period of Xq = 4 cm and an assumed undulator parameter K" = 1, (10) yields 
an electron beam energy 7 mc ^=100 MeV. The duty factor for this example is 
D = (3 ps)(500 MHz) = 0.0015. The average current \sT= Dt = 0.9 A which is much 
higher than typical FELs. To calculate the laser power which must leave the ship’s 
FEL to destroy the missile at a range of z = 8 km, the following steps are taken. The 
power required to destroy the missile is = 226 kW. The relationship between the 
power leaving the ship and the power arriving at the missile is = Pship^~^- 

Therefore, the laser power which must leave the ship is Pghip = 338 kW. The electron 
beam peak power is P = lV = (600 A )(100 MeV) = 60 GW which gives an average 
electron beam power of P =TV = 90 MW. The efficiency necessary to destroy an 
ASM is q = PghiplP = 0.004. An approximation for the natural single-pass efficiency of 
an FEL is q* = 1/(2A/) = 0.02 which is better than the required efficiency. Therefore, 
this example gives the required power density to destroy a missile. This type of 
calculation is the main driving factor in determining the FEL parameters for a possible 
weapons class FEL. 
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B. REVIEW OF MW WORKSHOP 


One threat to all Navy ship’s is an attack with an anti-ship missile (ASM). 
These weapons are part of the arsenal of many countries around the world, so the 
threat is real. These missiles are very fast and have a small radar cross section which 
makes them hard to detect. The Navy is performing more operations where a sole 
ship is sent to do a mission. An example is current anti-drug operations where a 
single ship operates close to shore for weeks at a time. Each ship must be capable of 
defending itself from ASM’s. For these reasons, an ASM is not something a battle 
group can be expected to defend individual ships against. Today, almost every US 
Naval ship has one or more PHALANX guns on board, the purpose of which is to 
defend the ship against an ASM attack. This is an illustration of how real the ASM 
threat has become. 

The current close-in-weapon-system (CIWS) weapon system has been tested 
against actual missiles. The author’s opinion is the CIWS is very capable of 
destroying an ASM. However, there is still a major problem with the PHALANX gun. 
The ASM is destroyed at such a close range that most of the debris, and a significant 
amount of energy, strikes the ship [Ref. 14]. While this may or may not be better than 
the intact ASM striking the ship, the ship is ultimately damaged by the ASM, which is 
the intent of the enemy. 

In an effort to replace this current weapon system, the Navy is exploring the 
use of an FEL as a ship self-defense (SSD) weapon. The range at which an ASM can 
be engaged increases from a few hundred meters to about 5-10 kilometers when a 
laser is used in place of the PHALANX gun. It has been estimated that it would take 
about 3 seconds for a laser beam to burn a hole in a missile cone. Because of the 
missile’s high speed, this hole would cause the missile to self destruct as a result of 
the reduction of aerodynamic properties and the ensuing forces this causes. Even for 
a missile traveling at Mach 3, if it is engaged at 5 km it will only have closed to a 
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range of approximately 2 km after the 3 second duration. Consequently, an FEL 
would keep the ship from being damaged by the ASM, thereby defeating the enemy’s 
attempt. 

An FEL is of interest for many reasons. A high power laser, on the order of a 
MW, will be required to destroy an incoming ASM. It is believed by commercial 
industry and the scientific community that high average powers can be achieved with 
FELs. Any SSD laser system will have to propagate in a maritime atmosphere where 
there are certain windows of propagation available. Being able to design the laser 
weapon system to a very specific wavelength is advantageous. The optical 
wavelength X of an FEL can be chosen by (10). If actual tests show that the desired 
wavelength is slightly different than originally believed, then modifying an FEL to adapt 
to the new wavelength is a relatively minor task. This ease of re-design is true for 
FELs but not for other lasers. For these reasons, the Navy is considering an FEL for 
ship self-defense. 

Currently, an infrared kilowatt FEL is under construction at Thomas Jefferson 
National Accelerator Facility (TJNAF) in Newport News, VA. This project, which has 
military as well as industrial applications, seeks to prove that a high average power 
FEL is a real possibility. In April 1996, the KW FEL project was started with US Navy 
funding. The Navy believes this laser will provide the Navy and the scientific 
community enough knowledge to be able to advance the FEL technology into a 
compact weapon design for ship self-defense. 

The Navy sponsored a MW class FEL concepts workshop at TJNAF on 
September 24-26, 1996. There were 57 attendees from government, national labs, 
universities, and industry. The tasking statement for this workshop was, "The 
workshop is to develop material for incorporation in a SSD FEL Systems Engineering 
Management Plan, and is to define baseline concepts and technology options for 
subsequent system performance and cost trades." [Ref. 15] While the Navy is not yet 
ready to construct a weapons class FEL, the meeting was designed to discuss the 
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feasibility of this idea, and to bring out some of the possible pitfalls which must be 
overcome. Out of the workshop came two possible designs which are discussed in 
the next section. 

C. WORKSHOP RESULTS 

Two conceptual designs came from the workshop. One design is an oscillator 
and the other is an amplifier. While these are very different approaches to a weapons 
class PEL, both are proposed to achieve the same goal. The goal is to produce 
enough average power to propagate energy through the atmosphere and burn a hole 
in a missile 5 to 10 km away. These two basic PEL designs evolved from a 
culmination of ideas presented at the workshop. 

The typical PEL oscillator design uses mirrors and small extraction efficiency to 
obtain the desired optical field strength. A simple schematic of the oscillator design is 
shown in Pigure 8. The PEL depicted in Pigure 8 has an overall length of 16.5 meters 



Pigure 8. This is a conceptual MW class oscillator PEL. 
Prom [Ref. 15]. 
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and a width of 5 meters. A uniform undulator only 1 meter long is responsible for 1% 
extraction efficiency. One unique characteristic is the fact that this design makes use 
of energy recovery. This concept involves extracting energy from the ymc^ = 100 MeV 
electron beam in the undulator. Then, before sending the electrons to the beam 
dump, they are decelerated. This deceleration allows the energy, which would 
otherwise be dumped, to be partially recovered. In the process, the electrons travel 
around a race track vacuum chamber. The same accelerator is used to decelerate 
and recover most of the left over beam energy. The design in Figure 8 has an 
average current 7 = 1 A. The input power is (1 A)(100 MeV) = 100 MW but 95 MW 
are recovered. With the output of 1 MW optical power this system requires a total 
input RF power of 6 MW. 

In the oscillator, the optical field increases gradually as the light bounces 
between the cavity mirrors. The optical field reaches high power which requires the 
mirrors to withstand high power density on their surface. If the power density on the 
mirrors is too high, mirror damage will occur. The power density placed on the mirrors 
can be reduced by designing the oscillator FEL with a short Rayleigh length [Ref. 16]. 
Depending on the design, the mirrors may or may not be at risk. 

An advantage of the oscillator design is that it appears to result in a more 
compact overall design. Space is critical on a ship so a smaller, but just as capable 
weapon, is a benefit. However, the overall volume of the oscillator and amplifier may 
be about the same. 

The fact that the oscillator design relies on energy recovery brings about both 
advantages and disadvantages. Extracting energy from the electrons leads to a 
smaller beam dump. On a ship, disposing of high energy electrons can be a problem. 
Radiation is given off during deceleration. This radiation can be harmful to both 
personnel and equipment if the beam is above 15 MeV energy. Shielding people and 
equipment from this radiation involves adding weight to a ship which is a great 
concern. So, the lower the energy of the electron beam entering the beam dump, the 
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better. To date, energy recovery has never been performed. One of the goals of the 
infrared kilowatt FEL at TJNAF is to demonstrate energy recovery. Another potential 
problem which stems from energy recovery is having to transport the high energy 
electron beam around a bend with the energy spread imposed by the FEL interaction. 
Electron beam transport is the subject of Chapter V. 

The regenerative amplifier design uses a longer undulator and relies on higher 
extraction efficiency to achieve high average power. A simple schematic of the 

s 

regenerative amplifier design is shown in Figure 9. The FEL depicted in Figure 9 has 



Figure 9. This is a conceptual MW class amplifier FEL. 
From [Ref. 15]. 


an overall length of 23 meters. A 6 meter long undulator is used to achieve 20% 
extraction efficiency from the 100 MeV electron beam. The undulator is uniform during 
the first half and tapered during the second half in order to achieve the higher 
extraction efficiency. The regenerative amplifier uses mirrors to feedback a small 
amount of the optical beam power which is then focused at the beginning of the 
undulator. At the end of the system, the high energy electron beam is diverted to a 
beam dump. The design in Figure 9 has an average current T = 0.05 A. This requires 
approximately 5 MW of input RF power and provides 1 MW of optical power in the 
infrared (IR) just as the oscillator design. 
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The amplifier design has advantages and disadvantages. This design does not 
use energy recovery so there is no requirement to bend the high energy beam. The 
bending may cause deterioration of the electron beam quality. Without energy 
recovery the electron beam is dumped at higher energy, which is a disadvantage. 
Even though there is no energy recovery in this design, the required RF power is 
similar to that in the oscillator design. Higher extraction efficiency requires a high 
quality electron beam and producing a high quality electron beam can be difficult. 

The two proposed MW designs include approximate parameters which describe 
them. These parameters are being improved based on knowledge obtained through 
computer simulations and experiments. The parameters which resulted from the MW 
workshop for both the oscillator and amplifier designs are listed in Table 1. The two 


Parameters 

Oscillator 

Amplifier 

Beam Energy 

100 MeV 

100 MeV 

Average Current 

0.9 A 

0.05 A 

Extraction Efficiency 

1% 

20% 

Undulator Type 

Linear 

Tapered 

Undulator Parameter 

1.0 

1.0 

Undulator Length 

1 m 

6 m 

RF Power 

6 MW 

5 MW 

Optical Power (IR) 

1 MW 

1 MW 


Table 1. These are the MW ship self-defense parameters 
for the two proposed FEL designs. 


designs have the same beam energy and they propagate the same wavelength 
radiation. They are different in most all other parameters. Typical FEL parameters are 
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listed throughout the theory chapter, Chapter 2. When comparing typical parameters 
with those for the MW designs in Table 1, the only real difference is found in the 
average current. The MW designs have a much higher average current, hence the 
resulting higher average power. Otherwise, the parameters listed are similar to any 
other FEL in the world. 
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IV. MW REGENERATIVE AMPLIFIER FEL 


A. DESIGN PARAMETERS 

In the previous chapter, two proposed MW FEL designs are discussed and the 
parameters for the two designs are listed in Table 1. These parameters resulted from 
the Navy’s MW meeting held in September 1996. Improvements have been made to 
each of the designs. Much of the work to optimize the oscillator design was 
conducted by the FEL group at NFS through computer simulations, the results of 
which are in Ref. 17. The regenerative amplifier design has been proposed by LANL 
and Boeing with the analysis being performed at NFS. An analysis of the regenerative 
amplifier parameters is the topic of this chapter. 

The design parameters for the MW regenerative amplifier (Figure 9) have 
evolved into two cases [Ref. 18]. Each is intended to provide the average optical 
power required to destroy an incoming anti-ship missile. The following parameters are 
common to both cases: the electron beam energy of 100 MeV, the normalized 
emittance of 6jr mm-mrad, and the rms energy spread which is Ay/y = 0.02%. In 
addition the electron beam energy (chosen to reach the 1 pm optical wavelength), the 
electron beam radius (0.17 mm), peak current (400 A), and pulse duration (20 ps) are 
common features. In each case, 4 m and 6 m long undulators are used in which the 
first half is uniform and the second half is tapered. These regenerative amplifiers will 
feedback on the order of 0.1% to 1% of the optical power and require an extraction 
efficiency of approximately 10% to 15%. These parameters are common to the two 
cases. 

In case 1, the undulator parameter K = 1.49, the undulator period ?io = 2.5 cm, 
and the gap is 5.3 mm. The gap is the distance between the upper and lower portion 
of the undulator and affects the strength of the magnetic field at the midpoint of the 
gap. The optical beam radius in the uniform section is 0.3 mm and there are /V = 160 
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undulator periods. The single pass gain is expected to be about 1000. Case 2 is 
more aggressive because it has a smaller gap size of 2.3 mm, a larger undulator 
parameter /<"= 1.71, and a shorter undulator period Xq = 2.0 cm. The optical beam 
radius in the uniform section is 0.25 mm and there are N = 200 undulator periods. 
This more aggressive design predicts single pass gain around 14,000. These 
parameters are listed in Table 2. 


Parameters 

Case 1 

Case 2 

Beam energy 

100 MeV 

100 MeV 

Beam radius 

0.17 mm 

0.17 mm 

Pulse duration 

20 ps 

20 ps 

Peak current 

400 A 

400 A 

K 

1.49 

1.71 

Xq 

2.5 cm 

2.0 cm 

N 

160 

200 

X 

1.045 pm 

1.045 pm 

Extraction efficiency 

> 15% 

> 15% 

Single pass gain 

~ 1000 

~ 14,000 


Table 2. Parameters for the MW regenerative amplifier designs 
developed at LANL. 


B. SIMULATIONS AND RESULTS 

Computer simulations are used as an evaluation tool. One-dimensional 
simulations are used to describe a single pass of the optical wave through the 
undulator. The regenerative amplifier feeds back from 1% to 0.1% of the optical 
power, so we require the simulation achieve gains of G = 100 to G = 1000 in order to 
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operate in steady-state. The initial optical field aQ is varied to achieve the desired 
gain. These simulations attempt to represent steady-state operation of the 
regenerative amplifier. The amount of taper in the undulator is determined by 5, the 
dimensionless artificial acceleration defined just above (36). The tapered undulator 
causes a shift in the phase velocity v away from resonance for the untrapped 
electrons. The electrons inside the separatrix are trapped and therefore responsible 
for gain. The taper turns on at the most efficient point along the undulator in order to 
maximize extraction efficiency. These two parameters ao and 5 are varied until the 
maximum possible extraction efficiency is obtained for gains of 100 and 10^. The 
results for each case are given below. 

1. Case 1 

The parameters for case 1, listed in Table 2, are used to derive dimensionless 
quantities. These dimensionless parameters are input into the computer simulation 
and the values of Bq and 5 are varied to optimize efficiency. The dimensionless 
current density is / = 1500 with A/ = 160 undulator periods. The electrons start on 
resonance, Vq = 0, to maximize high gain consistent with Figure 5. The small 
Gaussian spread oq = 0.4 represents an energy spread of 0.02% as seen from (34). 
In this simulation, the taper 5 = 70n turns on at t = 0.5. A large optical phase shift 
A<j) = 2n, typical of a tapered undulator with high gain and strong fields, is similar to 
that of Figure 7. The best result for case 1 is final gain G = 1100 and extraction 
efficiency q =4% given by ao = 15 and 5 = 70k shown in Figure 10. Increasing Bq 
tends to cause a decrease in gain while decreasing Bq tends to increase gain. 
Increasing or decreasing 6 primarily modifies the efficiency. Increasing 5 beyond the 
optimum value tends to trap fewer electrons and reduce efficiency. Decreasing 5 
lower than the optimum value does not make full use of tapering and also reduces 
efficiency. Gain is defined as 
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(37) 


G = 


a 


2 

0 


where la^ is the final field and Bq is the initial optical field strength. Assuming large 
gain G » 1 gives 

df = ao^ . (38) 


In case 1, a^ = 500 which means that if approximately 3% of the final optical field 
(0.1% of the power) is fed back and focussed at the beginning of the undulator, this 
system could operate in steady-state. While this feedback is reasonable, the 
extraction efficiency is substantially lower than the desired 15% indicating potential 
problems with this design. 


*** fEL Phase Space Evolution *** 
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Figure 10. Phase space evolution for the regenerative 
amplifier, case 1, with a 4 meter undulator. 
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In an effort to obtain ti = 15%, it was mutually decided by LANL and the FEL 
group at NFS to increase the length of the undulator in this design from L = 4 m to 
L = 6 m. The dimensionless parameters are again calculated to be used in the 
simulation. The increased undulator length is responsible for the increased values of 
current density to / = 5000 with N = 240 undulator periods and energy spread 
0(3 = 0.6. High gain j » n, dictates vq = 0 once again. In this simulation, the taper 
turns on at T = 0.3 and the stronger taper causes a larger difference in phase velocity 
V between electrons which are trapped and those which are not. This difference is 
seen when comparing the v versus ^ portion of Figures 10 and 11. The optical phase 
shift also occurs earlier in Figure 11 than the phase shift of Figure 10. The same 
previous steps are taken and the best results are final gain G - 1000 and extraction 
efficiency reached t] = 8.6%. This result is represented in Figure 11 where ao = 51 
and 5 = 150 tc. The maximum efficiency for any value of gain in Figure 11 is -n = 8.6%. 
As stated previously, increasing Bq tends to cause a decrease in gain while 
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Figure 11. Phase space evolution for the regenerative 
amplifier, case 1, with a 6 meter undulator. 
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decreasing causes gain to rise. When 6 is altered from its optimum value fewer 
electrons are trapped decreasing efficiency. The values of gain and initial optical field 
strength give Sf ~ 1630 so the required feedback is still approximately 3%. Even 
though gain is at the desired level, the extraction efficiency r\ is still less than the 
predicted 15%. 

In another effort to reach 15% efficiency, the final gain is decreased from 1000 
to 100. This is done by increasing both and 6. Increasing ao represents more 
feedback and stronger fields throughout the undulator, which leads to lower gain. It 
was expected that increasing a^ would increase the trapping efficiency and therefore, 
the overall efficiency. The optimum case yielding the greatest efficiency is shown in 
Figure 12 where ao = 122 and 5 = 122jt. The other parameters are the same as those 
in Figure 11. The taper starts at i: = 0.23 which is similar to the taper start in the 


*** FEL Phase Space Evolution *** 
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Figure 12. Phase space evolution for the regenerative 
amplifier, case 1, with a 6 meter undulator and gain 
G = 100. 
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previous figure. The result is G = 163 and ri = 7.8%, actually lower than the efficiency 
obtained when G > 1000. Increasing or decreasing either Bq or 5 has the same 
results as in the previous two trials. The reason for the decrease in efficiency here is 
that fewer electrons are trapped when the gain is lower. So, for case 1, the optimum 
results are represented by Figure 11 where ao = 51, 5= 1507c, 6 = 1025 and 

Ti = 8.6%. 


2. Case 2 

The same procedures followed in case 1 are now followed for case 2. The 
parameters used to develop dimensionless quantities are taken from the beginning of 
this chapter and Table 2. The goals for this design are large gain of G > 14,000 and 
high efficiency of ti = 15%. In the simulation, the dimensionless current density 
j = 3100 with N = 200 undulator periods, and energy spread <Sq = 0.5. As in case 1, 
the maximum efficiency found in computer simulations is lower than predicted. The 
optimum values found are given in Figure 13 where ao = 7 and 5=1107c. These 
results are 6 = 15,400 and t) = 4.7%. As in case 1, increasing Bq tends to cause gain 
to decrease and decreasing Bq results in increased gain. Changing 6 from the 
optimum value 110k results in a decrease in efficiency. When 5 is too high, fewer 
electrons are trapped reducing efficiency. When 6 is too low, the taper is not being 
taken full advantage of which also reduces efficiency. The taper 5 = 110k starts at 
X = 0.5. The calculated value for the final optical field Bf = 870 corresponds to about 
1% of the optical field being fed back to the beginning of the undulator. 

To obtain maximum efficiency, the goal of G > 14,000 is dropped. Varying 
both Bq and 5, the maximum efficiency is found at ao = 30 and 5 = 1 10k. The result is 
given in Figure 14 where gain is G = 1177 and efficiency is q = 6.6%. The taper in 
this simulation starts at x = 0.36. The increase in ao from 7 to 30 decreases gain so 
that the final field is Bf - 1029. If 3% of the field is fed back, this design could run in 
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Figure 13. Phase space evolution for the regenerative 
amplifier, case 2, with a 4 meter undulator. 


steady-state. However, this FEL does not produce the predicted extraction efficiency 
of 15%. 

As in case 1, the undulator length is increased to 6 meters to try and obtain a 
higher efficiency. The longer undulator increases j to 10,400, increases the number of 
undulator periods to A/ = 300, and increases Cq to 0.8. A stronger taper 5 = 2457: is 
used to obtain the maximum efficiency. The taper starts earlier at x = 0.23. The 
stronger taper causes a larger difference in phase velocity v between the trapped and 
untrapped electrons. The gain is also changed to G ^ 1000. The result of these 
simulations is shown in Figure 15 where G = 1554 and ti = 12.2%. The increased 
undulator length gives a higher value of the final optical field Bf = 3150. Therefore, 
with the same feedback of 3% of the field as in Figure 14, the value for the initial field 
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Figure 14. Phase space evolution for the regenerative 
amplifier, case 2, with a 4 meter undulator and no 
particular required gain. 


is Sq = 80. As with the previous two trials in this case, altering the optimum values of 
ao and 5 has the same results on gain and efficiency as before. This case nearly 
achieves the required goal of 15% efficiency and is the most promising. The reason 
for higher efficiency here is that a higher percentage of electrons become trapped with 
G = 1000 

As with case 1, the required gain is reduced from 1000 to 100 to try for higher 
efficiency. Again, the highest efficiency is found with G > 1000 instead of G ^ 100. 
The same parameters are used as in Figure 15 except for Sq and 5. Also, the taper 
starts at X = 0.17 which is a little earlier than in the previous figure. The final gain is 
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G = 139 and the efficiency is ti = 10% with initial field Sq = 243 and taper 5 = 17571, as 
shown in Figure 16. 
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Figure 15. Phase space evolution for the regenerative 
amplifier, case 2, with a 6 meter undulator and required 
gain of 1000. 
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★★★ FEL Phase Space Evolution *♦* 

j=10400 a =243 V =0 N=300 

o o 

0^=0.8 6=1757C T1=10% G=139 



-n;/2 c 37C/2 Ox 1 


Figure 16. Phase space evolution for the regenerative 
amplifier, case 2, with a 6 meter undulator and required 
gain of 100. 


C. FUTURE ANALYSIS 

One-dimensional computer simulations indicate that the two cases for the MW 
amplifier design will not yield the predicted extraction efficiency of t) = 15%. In the 
previous section, variations of the two cases gave values for the efficiency ranging 
from Ti = 5% to q = 12%. It remains to be determined if 15% efficiency is actually 
necessary or if a somewhat smaller value is acceptable in the FEL weapon design. A 
lower efficiency would require higher beam current to achieve the same optical power. 
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The resulting values for gain and extraction efficiency obtained in the previous 
section may actually be over estimated, as they do not include the effect of diffraction 
of the optical beam. The definition for the filling factor given in Chapter II is 
F = r^/WQ . The actual values used were F = 0.3 for case 1 and F = 0.4 for case 2, 
obtained from the original definition of the filling factor. If diffraction of the optical field 
occurs, then this definition of the filling factor is incorrect. When diffraction is present, 
the optical field will expand in area and decouple from the electron beam. From 
Chapter II, the dimensionless Rayleigh length is Zq = kwq/{LX). Case 1 gives a value 
of Zq = 0.07 where the optical beam radius in the uniform section is 0.3 mm, the 
undulator length L = 4 meters, and the optical wavelength X = 1.045 pm. For case 2 
Zq = 0.05. With such a short Rayleigh length, the optical beam would naturally expand 
in radius by a factor of = 20 by the end of the undulator. The optical wavefront would 
effectively uncouple from the electron beam 10% to 20% along the undulator, greatly 
reducing gain and efficiency. This FEL design relies heavily on optical guiding, 
discussed at the end of Chapter II, to keep the optical beam focussed near the 
electron beam. Therefore, further analysis of diffraction is suggested. 
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V. ELECTRON BEAM TRANSPORT 


A. REVIEW OF TRANSPORT WORKSHOP 

In chapter III, the results of the Navy’s ship self defense FEL workshop are 
discussed and two possible MW FEL designs are compared. The two designs are the 
oscillator design and the regenerative amplifier design which both have strengths and 
weaknesses. One characteristic specific to the oscillator design is energy recovery. 
While energy recovery is an advantage, it requires transporting a high energy bunched 
electron beam around four 180° bends. Imposing this bend on such a beam may 
cause coherent synchrotron radiation (CSR) to decrease beam quality. This effect has 
sparked Navy interest in CSR. 

The Navy is not the only party interested in these results. Two other 
communities are commercial industry and high energy physics. A high-average-power 
FEL is of interest to commercial industry for use in surface processing of materials with 
ultra-violet (UV) light. This application is believed to have global markets totaling 
hundreds of billions of dollars [Ref. 19]. Energy recovery is a must for industrial 
applications because it drives down the cost of light produced by increasing overall 
system efficiency. It also eliminates the massive radiation management program 
required by a high energy beam dump. The application in the field of high energy 
physics is in particle beam colliders. Therefore, there are many reasons for gaining 
understanding of GSR’s effect on beam quality. 

The anticipated problem is that of a short electron bunch experiencing an 
acceleration and radiating coherently. The wavelength of the radiation is slightly 
greater than the bunch length, typically 1 mm to 1 cm. This radiation is termed 
"Coherent Synchrotron Radiation." The problem with CSR is that it causes an energy 
spread in the longitudinal direction and an emittance growth in the transverse direction 
which results in a degraded electron bunch. In a degraded bunch, the change in 


45 


energy causes electron trajectories to change. If severe enough, some high energy 
electrons could actually hit the vacuum cavity walls. Also, a degraded bunch could 
make energy recovery difficult since an energy spread decreases the effectiveness of 
beam transport in the recovery process. Although there are currently several theories, 
the physical processes involved in bending the beam are complex and not well 
understood. While the prediction of CSR dates back to around 1950, the effects of 
CSR on emittance have never been measured. 

As a result of the Navy’s interest in this phenomena, a prototype beam-bending 
problem was drafted by the FEL group at NFS. The problem included a micropulse 
common to most of the proposed MW designs presented at the TJNAF meeting. The 
prototype micropulses’ peak current is 600 A with a full width, half maximum pulse 
duration of 3 picoseconds, and an electron beam energy of 100 MeV. The total bunch 
charge is 2.0 nC with a Gaussian longitudinal distribution. The nominal vacuum 
chamber width is 5 cm. The emittance and energy spread entering the bend are 
assumed perfect with zero spread. There is no specific bend design other than a 
continuous bend and one that is approximately achromatic and isochronous. What 
happens to this beam during a 180° bend with a bend radius of 2 meters? 

This problem was sent to approximately 20 scientists who have some 
experience with CSR for possible analysis. Next, an electron beam transport 
workshop was held at NFS on March 6, 1997. The purpose of this workshop was to 
examine the analyses and discuss the possible effects of transporting a high-energy 
electron beam around a 180° bend. Along with the analyses, the comments made by 
the workshop speakers were discussed. The organizations represented during this 
meeting were the Navy’s Directed Energy Office (SFAWAR), Los Alamos National 
Laboratory (LANL), Stanford Linear Accelerator Center (SLAC), Northrup Grumman, 
Thomas Jefferson National Accelerator Facility (TJNAF), and the Naval Fostgraduate 
School (NFS). 
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B. WORKSHOP CONCLUSIONS 


As expected, there was not total agreement among the analyses performed by 
the workshop participants. This stems from the many different CSR theories. 
However, a majority of the analyses gave similar CSR induced emittance growth. 
Overall, it is estimated that the CSR induced normalized rms emittance growth may be 
on the order of Ae^ = lOOn’s mm-mrad. The goal for a typical infrared PEL is to be 
around £„ = 2O71 mm-mrad, so a value of approximately IOO71 mm-mrad would result in 
a seriously degraded electron beam, and degraded PEL performance. 

Since the effects of CSR on emittance have never been observed, a CSR 
experiment would help benchmark theories. Tests are planned at LANL in the 
summer of 1997, and at TJNAP in January of 1998. As a result, the PEL community is 
on the verge of being able to compare theoretical models with measured results of 
experiments. One key in these observations is the actual accuracy of the emittance 
growth measurements. Measuring emittance growth accurately is a difficult task. 
Therefore, small emittance growths can sometimes go undetected. Por this reason, 
the participants of the workshop believe that there needs to be about a factor of two 
emittance growth in a system for actual detection to be feasible. In addition to the 
planned tests, workshop participants made a list of these and other possible 
experiments. Table 3 and Table 4 list the parameters for these possible CSR 
experiments. A simple analysis of these possible experiments is presented in the next 
section. As a result of this beam transport workshop, the Navy’s Directed Energy 
Office has stated an interest in finding other organizations which would be able and 
willing to aid the Navy in funding a CSR experiment. 

The parameters in Table 3 were obtained from each of the locations listed. 
There are two machines located at Los Alamos, one of which is an PEL. Both of 
these systems currently exist and are operational. Both the Princeton and Duke PELs 
are also operational. The TJNAP system is currently under construction and there are 
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Location 

Beam 

Energy 

(MeV) 

Peak 

Current 

(A) 

Average 

Current 

Normalized 

Emittance 

(71 mm-mrad) 

Energy 

Spread 

Pulse 

Length 

(PS) 

Beam 

Radius 

(fwhm) 

LANL (AFEL) 

17 

300 

0.5 A* 

<8 

<0.5% 

6-15 

0.5 mm 

LANL (Compressor) 

8 

1500 

10 nA 

10 

10% 

0.7 

200 p.m-5 mm 

Princeton 

9-14 

100 

0.22 A* 

10 

0.2% 

5 

2 mm-5 mm 

TJNAF (Chicane) 

42 

50 

5 mA 

13 

0.5% 

2.3 

0.4 mm 

(Recirculation Bend) 

42 

18 

5 mA 

14 

0.5% 

7 

2.5 mm 

Boeing 

110 

300 

100 mA 

15 

0.3% 

9 

0.5 mm 

Duke 

44 

40 

0.18 A* 

8-10 

0.3% 

2 

0.2 mm 

SLAC (Small Bend) 

14,350 

3400 

95 nA 

1.5 

0.02% 

0.23 

0.03 mm 

(Large Bend) 

14,350 

3400 

95nA 

1.5 

0.1% 

0.23 

0.03 mm 


* Average Current during macro-pulse. 

Table 3. Electron beam parameters for possible CSR experiments. 


two places in which emittance measurements will be made. The first is before and 
after a chicane, which is a bend used for bunch compression. The other place is 
before and after a 180° recirculation bend. There are two recirculation bends which 
allow this system to use the energy recovery technique. When looking for emittance 
growth in the recirculation bend, the undulator will be removed from the system so a 
strict CSR-induced emittance growth experiment, without undulator interaction, will be 
performed. Assembly of the Boeing PEL was stopped prior to commissioning trials 
due to lack of funding. The system at SLAC is still in the planning stages and is 
proposed to have two possible emittance growth measurement locations, each with 
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slightly different parameters. The beam radius is given as a full width, half maximum 
(fwhm) value. 


Location 

Undulator 

Periods 

N 

Undulator 

Wavelength 

^0 

Undulator 

Parameter 

K 

Bend 

Radius/Arc 

LANL (AFEL) 

100 

2 cm 

1.0 

- 

LANL (Compressor) 

- 

- 

- 

0.25 m/45° (a) 

Princeton 

73 

1.4 cm 

0.2 

0.2 m/90° (b) 

TJNAF (Chicane) 

40 

2.7 cm 

0.7 

1 m/22° (c) 

(Recirculation Bend) 

- 

- 

- 

1 m/180° (d) 

Boeing 

220 

2.18 cm 

1.31 

42 cm/180° (e) 

Duke 

47 

2.3 cm 

1.0 

43 cm/45° (f) 

SLAC (Small Bend) 

3330 

3 cm 

3.7 

25 m/13.6° (g) 

(Large Bend) 

3330 

3 cm 

3.7 

75 m/5.2° (b) 


Table 4. Undulator and beam transport parameters for possible CSR 
experiments, (a) Chicane, achromatic, (b) Two 45° dipoles separated by two 
quads, forms 90° achromatic bend, (c) Four 22° dipoles for bunch compression 
and another identical chicane for bunch compression following the wiggler. (d) 
Two separate 180° recirculation bends, achromatic and isochronous, (e) Four 
45° bends form a 180° bend, doubly achromatic and nearly isochronous, (f) 
Two 4-dipole chicanes, 11.25° per dipole, achromatic, (fir) Four bend chicane, 
3.4° per bend, achromatic, (b) Four bend chicane, 1.3° per bend, achromatic. 


Table 4 lists parameters for the undulators and the bends of the possible CSR- 
induced emittance growth experiments. The FEL at Los Alamos (AFEL) is a straight 
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system which does not bend the electron beam before the undulator. However, as an 
FEL, the undulator does bend the electron beam repeatedly. It is not known whether 
this back and forth bending causes cancellation of the CSR effect or not. Therefore, 
one of the possible experiments is to measure any CSR wavelength radiation 
produced by the undulator. The compressor at Los Alamos is not an FEL, therefore 
there is no undulator. The TJNAF recirculation bend experiment, as previously stated, 
will be conducted without an undulator in the system. Therefore, no undulator 
parameters are listed for the recirculation bend experiment. 

C. ANALYSIS OF POSSIBLE EXPERIMENTS 

A large CSR-induced emittance growth in a CSR experiment makes the effect 
easier to measure. Predicting this growth with computer simulations for each of the 
possible experiments can be a rather time-consuming task. Instead, a simplified 
analysis which examines steady-state CSR fields is conducted which enables the 
different experiments to be compared. This relative comparison is used to anticipate 
the experiments which may prove to be the most useful in learning about the CSR 
phenomena. 

In this analysis, the final normalized emittance e is given by 

8 = Veo + (Ae)^ (39) 

where Eq is the initial normalized emittance of the beam at the bend entrance, and Ae 
is the change in normalized emittance due to CSR fields. The wavelength of the CSR 
emitted by the electron bunch is on the order of the bunch length itself. This coherent 
radiation wavelength is given by 

X > 2kc (40) 

where o is the rms bunch length [Ref. 20]. For example, a pulse length c = 0.5 mm 
yields radiation X>k mm. Now, depending on the size of the vacuum chamber in 
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which the beam is transported, some or all of the CSR may be shielded by the walls of 
the cavity. For CSR fields to be unshielded by the vacuum chamber, the minimum 
chamber dimension should be much greater than 


1/3 


w * = \r{kc)^ 


(41) 


where w' is the shielding dimension, and R is the bend radius. The normalized 
emittance growth At is estimated using 



(42) 


where is the transverse rms beam size averaged along the beam path. Ip is the 
peak current in the electron pulse, Ia = 17,000 A is the Alfven current, and 0 is the 
bend angle of the beam path. This estimate of Ae assumes the electron bunch is rigid 
during the bend and the longitudinal CSR force is the main cause of emittance growth. 
Also, the transverse CSR force, which may contribute to emittance growth, is ignored. 
And, only single bunch effects are examined which means no effects of one bunch on 
another are included. [Ref. 21] 

Table 5 and Table 6 lists the parameters required in this analysis. These 
parameters are constructed based on Tables 3 and 4. In Table 5, the AFEL at Los 
Alamos is not included since there are no transport bends in this system. The electron 
beam sizes are given as average values instead of a range of values and have been 
converted from fwhm to rms. 

Table 6 lists parameters used to evaluate the different experiments. The 
values in the shielding chamber dimension column are calculated using (41). w is the 
actual minimum dimension of the vacuum chamber during the bend of each system. 
The right section of the table shows the estimated CSR-induced normalized emittance 
growth calculated using (42) and the final emittance calculated using (39). The 
percentage increase in emittance caused by CSR is also given. The three columns in 
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a 


R 

e 

Location 

Peak 

Bunch 

Transverse 

Bend 

Bend 


Current 

Length 

rms beam 

radius 

angle 



(rms) 

size (avg) 



LANL (Compressor) 

1500 A 

0.09 mm 

2.6 mm 

0.25 m 

45° 

Princeton 

100 A 

0.64 mm 

3.5 mm 

0.2 m 

CO 

o 

o 

TJNAF (Chicane) 

50 A 

0.29 mm 

0.4 mm 

1.0 m 

o 

CO 

CO 

(Recirculation Bend) 

18 A 

0.89 mm 

2.5 mm 

1.0 m 

180° 

Boeing 

300 A 

1.15 mm 

0.5 mm 

0.42 m 

180° 

Duke 

40 A 

0.26 mm 

0.2 mm 

0.43 m 

45° 

SLAC (Small Bend) 

3400 A 

0.03 mm 

0.03 mm 

25 m 

13.6° 

(Large Bend) 

3400 A 

0.03 mm 

0.03 mm 

75 m 

5.2° 


Table 5. Parameters for analysis of possible CSR experiments. 


this right-hand section all assume that there is no CSR shielding by the vacuum 
chamber, w » w*, which is not true for Duke and SLAC. 

Comparing w and w’ indicates the experiments located at Duke and SLAC 
may shield some or all of the CSR as a result of their minimum vacuum pipe 
dimension since w is not much greater than w' as required. On the basis of 
percentage emittance growth LANL, Princeton, Boeing, and SLAC have at least a 
factor of two emittance growth. Looking at all the parameters in Table 6 and 
remembering which FELs could support a CSR experiment in the near term, LANL 
Compressor and the Princeton FEL seem to be the best candidates. 
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Location 

w* 

Shielding 

chamber 

dimension 

w 

Minimum 

chamber 

dimension 

£o 

Initial 

emittance 

(Ttmm-mrad) 

Ae 

Emittance 

growth 

(Ttmm-mrad) 

8 

Final 

emittance 

(Tcmm-mrad) 

(e/8o) - 1 

% 

Emittance 

growth 

LANL (Compressor) 

0.27 cm 

2.0 cm 

8 

180 

180.2 

2150% 

Princeton 

0.93 cm 

3.4 cm 

10 

18 

21 

110% 

TJNAF (Chicane) 

0.94 cm 

5.0 cm 

13 

2 

13.2 

12% 

(Recirculation Bend) 

1.98 cm 

5.0 cm 

14 

12 

18 

29% 

Boeing 

1.76 cm 

7.0 cm 

15 

28 

32 

113% 

Duke 

0.66 cm 

1.3 cm 

10 

0.3 

10 

0% 

SLAC (Small Bend) 

0.61 cm 

0.2 cm 

1.5 

4 

4.3 

187% 

(Large Bend) 

0.87 cm 

0.2 cm 

1.5 

0.7 

1.7 

13% 


Table 6. Resulting parameters for analysis of possible CSR experiments. 


D. PRINCETON PEL 

One of the most readily available systems for a CSR experiment is the 
Compact Infrared PEL (ClREEL) located at Princeton University. This PEL was 
developed by Northrop/Grumman Corporation and Princeton. The system is 
operational and encorporates a 90° bend prior to the undulator as part of the electron 
beam transport path. Parameters of this PEL are listed in Tables 3-6. Emittance 
could be measured at the entrance and exit to the bend with only slight modifications 
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to the system [Ref. 22]. The analysis in the previous section approximately predicts a 
factor of two emittance growth caused by CSR in the 90° bend making a CSR 
experiment feasible. Depending on the accuracy of those measurements, this could 
be the first exact measurement of the CSR effects on emittance. For these reasons, 
the CIRFEL seems to be a fruitful and possibly inexpensive tool for performing a CSR 
experiment. 

Because of other priorities and funding issues, the CSR experiment has not 
already been performed. However, in an effort to further understand the CIRFEL, 
computer simulations of the FEL have been performed by Northrop/Grumman and the 
FEL group at NFS. Ongoing comparisons of these simulations hope to prove helpful 
in determining the limits of the CSR effect. 

Since CIRFEL is operational, the gain of the system has been roughly 
measured. Comparison of the predicted and measured gain may give insight as to 
whether or not CSR induced emittance gro\wth has occurred. Gain and emittance are 
coupled through beam quality. Large emittance gives poor beam quality which may 
cause lower gain. To calculate dimensionless parameters for the CIRFEL, some 
additional information about the electron beam is used. Northrop/Grumman reports 
that the full width, half maximum electron beam size in the vertical direction is 2 mm. 
In the horizontal direction the size of the beam varies. Horizontally, the beam is 5 mm 
at the undulator entrance, decreases to 2 mm at the center of the undulator, and is 
reported to be 10 mm at the undulator exit. To establish the effect of beam size on 
the FEL gain, the nominal sizes of 2 mm, 3 mm, and 4 mm are examined. Computer 
simulations at NFS yield final gain of 67%, 28%, and 20%, respectively, for the 
different beam sizes. Northrop/Grumman estimates gain around 30% which most 
closely corresponds to the NFS simulation with a beam size of 3 mm. The NFS 
simulations use initial emittance of IOtt mm-mrad and emittance growth in the bend of 
Ae = 1871 mm-mrad so the emittance seen at the entrance to the undulator is e = 2l7t 
mm-mrad. These values are taken from Table 6. 
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In conclusion, the analysis in this chapter predicts a factor of two CSR induced 
emittance growth and gain of approximately 30% for the CIRFEL. Because the actual 
initial emittance of the CIRFEL has not been measured, actual CSR induced emittance 
growth is not predicted. Measurements of gain and of emittance before and after the 
bend will be an excellent benchmark for this analysis and is recommended. 
Performing this experiment will help the Navy decide if the MW oscillator FEL design is 
a real possibility or if energy recovery is risky due to CSR. 
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